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Abstract 

We study possible crossing of the phantom divide line in a DGP-inspired F(R, (f>) 
braneworld scenario where scalar field and curvature quintessence are treated in a 
unified framework. With some specific form of F(R, <p) and by adopting a suitable 
ansatz, we show that there are appropriate regions of the parameters space which 
account for late-time acceleration and admit crossing of the phantom divide line. 
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1 Introduction 



Based on several astronomical evidences, our universe is currently in a period of positively 
accelerated expansion [1]. It is possible to interpret this late-time acceleration based on yet 
unknown component called dark energy in literature ( see for instance [2] with a comprehen- 
sive list of references therein). Also, it has been shown that such an accelerated expansion 
could be the result of a modification to the Einstein-Hilbert action ( for a recent review 
see [3]). On the other hand, DGP braneworld scenario has the capability to interpret this 
late-time acceleration via leakage of gravity to extra dimension in its self-accelerating branch 
[4]. For the first alternative, the simplest candidate for dark energy is the cosmological con- 
stant itself. However, it suffers from serious problems such as a huge amount of fine-tuning 
[2,5,6]. Within dark energy viewpoint, for a dark energy component, say X, its equation of 
state or equivalently equation of state parameter ui x = jj^-, determines both gravitational 
properties and evolution of the dark energy. Recent constraints on equation of state of dark 
energy indicate that px ~ —px and even that Ux < — 1- When we consider background 
evolution, there is no problem with ux < — 1. However, there will be apparent divergencies 
in perturbation when one crosses the phantom divide line, u>x — — 1- Accordingly, crossing 
of the phantom divide line, ux — — 1, provides a suitable basis to test alternative theories 
of gravity or higher dimensional models which can give rise to an effective phantom energy 
[6]. With these motivations, the issue of phantom divide line crossing has been investigated 
extensively in recent years ( see [6] and references therein). In principle, crossing of phan- 
tom divide line by the dark energy equation of state parameter at recent red-shifts, has two 
possible cosmological implications: either the dark energy consists of multiple components 
with at least one non-canonical phantom component or general relativity needs to be ex- 
tended to a more general theory on cosmological scales. Both of these conjectures have been 
studied extensively. For instance, curvature quintessence [7] is a fascinating proposal in this 
respect. Recently, phantom-like behavior in a brane-world model with curvature effects and 
also in dilatonic brane-world scenario with induced gravity have been reported [8]. Dark 
energy models with non-minimally coupled scalar field and other extensions of scalar-tensor 
theories have been studied widely some of which can be found in reference [9]. In the spirit 
of scalar-tensor dark energy models, our motivation here is to show that a general F(R, <p) 
DGP-inspired scenario can account for late-time acceleration and crossing of phantom di- 
vide line in some suitable domains of model parameters space. To show this feature, first we 
study cosmological dynamics of F(R, <fi) DGP-inspired scenario briefly. In the minimal case, 
motivated by modified theories of gravity, some authors have included a term of the type 
t§R n m the action [7,10]. This extension for some values of n has the capability to explain 



2 



late-time acceleration of the universe in a simple manner. It is then natural to extend this 
scenario to more general embedding of DGP inspired scenarios. The purpose of this paper is 
to perform such a generalization to study both late-time acceleration and possible crossing of 
phantom divide line in this setup. We focus mainly on the potential of the type V{4>) = Xcj) 11 
for non-minimally coupled scalar field. Our setup for minimal case predicts a power-law ac- 
celeration supporting observed late-time acceleration. In the non-minimal case, by a suitable 
choice of non-minimal coupling and scalar field potential, one obtains accelerated expansion 
in some specific regions of parameters space. While a single minimally coupled scalar field 
in four dimensions cannot reproduce a crossing of the phantom divide line for any scalar 
field potential [11], a non-minimally coupled scalar field account for such a crossing [6]. In 
DGP model, equation of state parameter of dark energy never crosses oo(z) = —1 line, and 
universe eventually turns out to be de Sitter phase. Nevertheless, in this setup if we include 
a single scalar field (ordinary or phantom) on the brane, we can show that equation of state 
parameter of dark energy can cross phantom divide line [12]. Crossing of phantom divide line 
with non-minimally coupled scalar field on the warped DGP brane has been studied recently 
[13]. Our purpose here is to obtain an extension of these non-minimal dark energy models 
within DGP-inspired F(R, <fi) braneworld scenario. We use a prime for differentiation with 
respect to R. An overdot marks differentiation with respect to the brane time coordinate. 

2 F(R, (j)) DGP-inspired Gravity 

We start with the following action 

3 2 i 

S=^- jd 5 (x) v ^^ + J d \x) v ^(^F(R, ( f>)--q^V, ( f ) VA-V( ( f>)+mlK + L m ), (1) 

where the first term shows the usual Einstein-Hilbert action in 5D bulk with 5D metric 
denoted by qab and Ricci scalar denoted by 3?. The second term on the right is a gener- 
alization of the Einstein-Hilbert action induced on the brane. This is an extension of the 
scalar-tensor theories in one side and a generalization of /(i?) -gravity on the other side. We 
call this model as F(R, <p) DGP-inspired scenario, y is the coordinate of the fifth dimension 
and we suppose that brane is located at y — . q^ u is induced metric on the brane which is 
connected to qab via = 5 IM A 5 U E 'qab- K is the trace of the mean extrinsic curvature of 
the brane defined as follows 

^ = \ £ !™: (\K»Ay=-e + [K,X=+) (2) 
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We denote matter field Lagrangian by L m (q^, ip) with the following energy-momentum ten- 
sor 



5q 



pv 



(3) 



The pure scalar field lagrangian is = ^q^V ^(pV u <f) — V((f>) which gives the following 
energy-momentum tensor 



= V^V^ - -<v(V0) 2 - q^V{<j)) 



(4) 



The field equations resulting from this action are given as follows 

JT^^B-\gAB^) +ml8 A ^ B v (R lll/ - l -q, v R)5{y) = 5 A »5 B \% v a- f lw + T^r)8(y). 

(5) 

In this relation T^ v = F ,mj,\ where T^ u is the energy-momentum tensor in matter frame and 



= p^f^- Also, Tf% rv is defined as follows 



7D ^ r l 

rpcurv _ 3_ 

" F'(R, 



-q^ ( F(R, 0) - RF'(R, 0) ) + ( F'(R, 0) ) ( q m q v(i - q, u q a p 



;af3 



(6) 



In the bulk, Tab = and therefore 



Gab = ^-ab — ~9ab^ = 



and on the brane we have 



1 T 

Guv Ruu ~zQu,vR 9 ' 

2 mi 



(7) 



(8) 



where r M „ = T^ u + T ^ + T™ rv . The corresponding junction conditions relating quantities 
on the brane are as follows 



limji^]^ = ^ 



e >+0 l 



ml 



y=0 



—^F'(R, 0) 

m% 



Rfiv — -QpuQ al3 Ra/3 



(9) 



A detailed study of weak field limit of this scenario within harmonic gauge on the longitudinal 
coordinates and using Green's method to find gravitational potential, leads us to a modified 
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(effective) cross-over distance in this set-up as follows ( see [14] for details of a similar 
argument) 



m 



Ud_F\ = fdF\ 



2ml\dRJ \dRJ 



DGP> 



(10) 



where 



DGP 



The gravitational potential in this scenario takes the following forms in 
two different extreme: 
For r <C £p 



U(r) 

and for r»f f 



(Mjp + M ) 
Q7im 2 3 F'(R,(j))r 



) qrt iy qn T* 

1 + ^-^-p + - F l ^ +0 ^p\ 



U{r) 



(Afy + M ) 



9f 2 f 4 



(12) 



Qit 2 m\r 2 

is effective mass of the scalar field which is non-minimally coupled to induced gravity, 
My, collectively shows the effective mass of other material fields minimally coupled to gravity 
and 7 = 0.577 is Euler's constant. The most important feature of this F(R, 0) DGP-inspired 
scenario is the fact that now cross-over scale is explicitly related to the induced curvature on 
the brane and non-minimal coupling of scalar field with this curvature. Since the dynamics 
of scalar field is described by V M 0V M — ^ + -tt^ = 0, one may explain this result 
as a spacetime variation of the Newton's constant. In this viewpoint, when F(R, <fi) varies 
from point to point on DGP brane, the crossover scale takes different values. Recent best-fit 
crossover scale is given by Idgp = [15] where Hq 1 ~ 3000Mj9c. Since ^ = 2m Jx F , 

if we choose lp ~ \-2QHq 1 it is possible to constraint this F(R, 0) scenario to be consistent 
with observational data. For instance, if we set F(R, 0) = a{(f>)R it is easy to show that a(0) 
can attain the value of ct(<f>) ~ 2.52 x IO^Hq 1 where we have assumed m 4 ~ 10 3 GeV and 
m 3 ~ lO^GeV. This argument has the potential to be used as a mechanism for obtaining 
reliable form of F(R, 0) in cosmo logical context. 



3 Cosmological Implications of the Model 

Embedding of FRW cosmology in DGP scenario (which deviates from general relativity at 
large distances) is possible in the sense that this model accounts for cosmological equations 
of motion at any distance scale on brane with any function of Ricci scalar. It is well-known 
that original DGP model explains late-time acceleration via leakage of gravity to extra 
dimension. On the other hand, in this model equation of state parameter of dark energy 
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never crosses the phantom divide line. Extension of DGP scenario to more generalization 
of the brane action may result in several new implications on cosmological ground. In this 
respect, it is interesting to know whether DGP-inspired F(R, <fi) gravity can account for 
late-time acceleration and especially phantom divide line crossing. With this motivation, 
in this section we investigate cosmological implications of our setup focusing firstly on the 
late-time acceleration. We assume the following line element 



ds 2 = q^dx^dx" + b 2 (y, t)dy 2 = -n 2 (y, t)dt 2 + a 2 (y, t)j ij dx l dx j + b 2 (y, t)dy 2 



(13) 



where 7^ is a maximally symmetric 3-dimensional metric defined as 7^ = 5^ + k where 
k = —1, 0, 1 parameterizes the spatial curvature and r 2 = XiX 1 . We assume that the scalar 
field depends only on the cosmic time of the brane. Choosing a Gaussian normal coordinate 
system so that b 2 (y,t) = 1, non-vanishing components of Einstein's tensor in the bulk plus 
junction conditions on the brane defined as 



■. — >+u 



2nm 



2 r 



n 



3m\ 



ml 



'dF\ 
^~dR> 



>dF 



X-4-- 



\im Q [d y a} y y 



(*) 



ml 



fdF\ / a 2 k\ 
\dR)\tfa~ + a) 



a 2 


ha k 


2n 2 a 2 


n 3 a 2a 2 


(tot)^ 








y= 


1 








VcLF\ p('°V 


- y=o 




XdRj 3ml . 



y=0 



y=o 



(14) 



(15) 



yield the following generalization of Friedmann equation for cosmological dynamics on the 
brane ( see [14] for machinery of calculations for a simple case) 



H 2 + 



k_ _ 1 

a 2 ~ 3m 2 3 F'(R, 0) 



P M + Po 



1 + e 



1 + 



Po 



ptot 



m§F'(i?,0)e o 



(16) 



where e = ±1 shows two different embedding of the brane, p — m ^F'(R <t>) anc ^ £ ° = 
3(^1 — a' 2 + kja 2 is a constant ( with a' = ^). Total energy density and pressure are defined 



as p^ = p + p + p c 



p(tot) _ p _|_ _|_ pcurv reS p ec tively. The ordinary matter on the 
brane has a perfect fluid form with energy density p and pressure p, while the energy density 
and pressure corresponding to non-minimally coupled scalar field and also those related to 
curvature are given as follows 



P v = 



1 , 



dF 



-<p 2 + n 2 V{4>)-Q—H<P 



(17) 
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also 



1 ^ _ y + _2_d*\ v _ VLj,) + 4^^0 + — — 2 
2n 2 n 2 <i0 n g?0 n 2 n 2 d(f) 2 



y=0 



ml (\ 



F'(R,</>)\2 



F(R,(f)) - RF'{R,<f)) -3RHF"(R, </>)), 



(18) 



(19) 



P™™ = [2RHF"(R, 0) + i?F"(i?, 0) + i? 2 F w (i?, 0) - - 0) - 0)] J .(20) 

where H = is the Hubble parameter. Ricci scalar on the brane is given by 



k 1 
R = 3^ + 



a 2 ri 2 



6- + 6(-) -6- 



a n 
n 



In this setup, non-minimal coupling of scalar field and induced gravity leads to no-conservation 
of effective energy density on the brane 



+ 3H(p^ + v ^)) = 6(^)0(# 



d(j>- 



+ 



k 



(21) 



It is easy to show that for a minimally coupled scalar field on the brane, this setup yields 
a late-time accelerating universe in a fascinating manner [10]. The evolution of the scalar 
field on the brane for a spatially flat FRW geometry, (k — 0), is described by the following 
equations 



1 



(tot) 



and 



H 3m 2 F'( J R,0) P 

dV{4>) dF{R,<j>) 



+ 3#0 + 



dip 



d(j) 



(22) 



(23) 



and finally 



d(0,t) 1 
a(0,t) " ~~ §m\F'{R,$) 



p tot + 3p 



lot 



(24) 



These equations show that essentially embedding of FRW cosmology in DGP setup is pos- 
sible. Now, after a brief study of gravitational and cosmological implications of our setup, 
we investigate late-time acceleration and crossing of the phantom divide line in this setup. 
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3.1 Accelerated Expansion 



Standard DGP model itself has the capability to explain the late-time acceleration of the 
universe via leakage of gravity to extra dimensions without any additional mechanism [4]. 
In our F(R, <fi) DGP inspired model, it is interesting to see whether there is any room for 
explanation of this late time acceleration. The viability of this question lies in the fact that 
generally with non- minimal coupling it is harder to achieve accelerated expansion [16]. We 
first obtain a necessary condition for the acceleration of the universe in F(R, 4>) DGP-inspired 
model. Then we use a simple ansatz to clarify our general equations. Suppose that scalar 
field is the only source of matter on the brane so that p = for other matter candidates. 
The necessary condition for acceleration of the universe is (p^ + p curv ) + 3(p^ + p c «™>) < o 
because in our F(R, 0) setup scalar field and curvature are correlated. Therefore we obtain 

Using the Klein-Gordon equation (23), this relation can be rewritten as follows 



This is a general condition for acceleration of the universe in our F(R, <f)) DGP inspired 
model. It is a complicated relation and one cannot achieve an explicit intuition of acceleration 
in this setup. For simple forms of F(R, <fi) we can find relatively simple conditions that can 
be explained more explicitly. As a simple example and following Faraoni [16], if we set 
F{R,4>) = |(1 — i4> 2 )R which gives conformal coupling of scalar field and Einstein gravity, 
equation (25) under assumption of weak energy condition p^ > 0, gives V-f(j)^ > 0, where 
we have assumed £ < | [16]. For instance, if we set V(4>) = A0 T/ , we find rj < Using the 
ansatz a(t) m At v and <f>(t) « Bt'^, equations (22) and (23) with F{R,<j>) = \{l - £4> 2 )R 
and p — 0, give £ > ^ with positive and real v and considering terms of order 0{t~^~ 2 ). 
In this case £ is restricted to the interval < £ < | and since equation (23) with these 
ansatz gives p? + (1 — 3u)p + 12£z/ 2 + 6£z/ = 0, positivity and reality of solutions for p gives 
(9 - 48£)^ 2 - (6 - 24£)z/ + 1 > 0. For £ = we find v = (4 ± v / ll)/5 which gives a power- 
law accelerated expansion for positive sign. Therefore, a suitable fine-tuning of non-minimal 
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coupling provides the possibility of late-time accelerated expansion. Based on a dark energy 
model, to have an accelerated universe, the value of conformal NMC should be restricted 
to the interval 0.146 < £ < 0.167 [17]. On the other hand, current experimental limits 
on the time variation of G constraint the nonminimal coupling as — 10~ 2 -< £ -< 10~ 2 [18]. 
Solar system experiments such as Shapiro time delay and deflection of light have constraint 
Brans-Dike parameter to be uj B d > 500 which leads to the result of |£| ^ 2.2 x 10~ 2 for 
non- minimal coupling [18]. However, other approaches lead to different constraints on the 
value of conformal coupling [19]. As a second and more general example with R n correction, 
we set F(R, <fi) = \(1 — ^4> 2 )^R n . For spatially flat FRW geometry the Ricci scalar is given 
by 

R = 6- + 6(-) 2 . (26) 

Using the above proposed ansatz for a(t) and 0(t), equations (22) and (24) can be rewritten 
as the following explicit time- dependent form 

3m 2 z/ 2 £ n(12z/ 2 - 6i/) n_1 (l - fB 2 r 2 ") = 

5yr 2 ^ +2n - 2 +2A5r,r w+2 ™-^ 

71 (27) 

-3m 2 3 v(is - l)£ n(12u 2 - Qu) n - l r 2n+2 [l - £B 2 r 2fI ) = 
25yr 2 ^-6£5 2 /i4(12z/ 2 -6z/) n r 2n - 2 ^ 

+Xr]Bn- m ) + 2\Bn~ m+2 + 6m 2 Jn - l)(2n - v - 1) - m 2 (- - l)(12z/ 2 - 6u). (28) 
' n 

Our aim is to see whether these equations account for positively accelerated expansion with 
suitable values of parameters. Since the parameter space of the problem is complicated, 
analytical solutions of these equations have no obvious interpretation. Therefore, we try to 
see how a suitable choice of parameters leads to viability of accelerated expansion in this 
setup. Our strategy is to see how with different choices of parameters in equations (27) 
and (28), equality in this equations is preserved. Firstly, we choose v > 1 in the favor of 
late-time positively accelerated expansion. We find appropriate values of other parameters 
in this parameters space as shown in table 1 and 2. We plot each side of equation (27) as a 
separate function of time for some values of parameters. Existence of intersection point for 
these functions shows the viability of values attributed to parameters to preserve equality. 
The same procedure for both sides of equation (28) is done too. In summary, to have 
accelerated expansion, we need v > 1, so choosing v > 1 we find some possible values of 
the other parameters to have at least one intersection point in the graph and therefore to 
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Table 1: Some choices of parameters for equation (27) to be preserved. The last column 
shows the coordinate of intersection points and equality of two sides of equation (27) with 
corresponding choices of parameters. 



V 


n 


V 




(a, b) 


2 


1.4 


1.2 


-0.166 


(1.4, 16.4), (7.6, 15.2) 


2 


1.4 


1.2 


0.166 


(7.6,15.2) 


4 


1.4 


1.2 


-0.166 


(1.4,16.5) 


4 


1.4 


1.2 


0.166 




2 


1.5 


2 


-0.166 


(1.6,114) 


2 


1.5 


2 


0.166 




4 


1.5 


2 


-0.166 


(1.6,114) 


4 


1.5 


2 


0.166 





be consistent with accelerated expansion. Secondly, choosing some specific values for other 
parameters, we obtain possible values of v. If there is no intersection between two graphs, 
it means that with that choice of parameters we will not arrive at v > 1 corresponding to 
accelerated expansion. If there is intersection but with v < 1, the solution will be decelerating 
ruled out by observational data. The results of these analysis are shown in table 1 and 2 
and also in figures 1, 2, 3 and 4. Note that these are only some specific examples and the 
actual parameter space is very large and complicate. However, the least advantage of this 
analysis is the fact that essentially late-time acceleration can be explained naturally within 
our setup. We see for instance that with rj — 4, n — 1.4, and £ = 0.166 there is no possibility 
to have accelerated expansion with v — 1.2 since there is no intersection point in the graph. 
On the other hand, accelerated expansion with v — 2 requires 77 = 2 or = 4, n = 1.5, 
and £ = —0.166 for instance. As a common role in this analysis, we have noticed that from 
equation (27), only for positive values of n with 1.37 < n < 2 and for positive or negative 
£ there are intersection points and therefore accelerated expansion. These analysis, though 
very simple and especial, show that F(R, 4>) DGP-inspired scenarios essentially account for 
positively accelerated expansion. 

Note that in minimal case where F(R, <fi) oc R n one finds from (22) and (24) [20] 

i/[i/(ra-2) + (n-l)(2n-l)] = (29) 
v[v{n — 2)+n 2 -n + l] =n(2ra-l)(n-l) (30) 
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Table 2: Some choices of parameters for equation (28) to be preserved. The last column 
shows the coordinate of intersection points and equality of two sides of equation (28) with 
corresponding choice of parameters. Note that these are only some example and the actual 
parameter space is very wide. 



V 


n 


V 




M) 


2 


-0.37 


1.005 


-0.166 


(63.4,46.4) 


2 


-0.37 


1.005 


0.166 


(122,281) 


4 


-0.37 


1.005 


-0.166 


(92.4, 128) 


4 


-0.37 


1.005 


0.166 


(92.4, 128) 


2 


-0.5 


1.2 


-0.166 


(16.9,52) 


2 


-0.5 


1.2 


0.166 


(20.6,96.3) 


4 


-0.5 


1.2 


-0.166 


(19,69.4) 


4 


-0.5 


1.2 


0.166 


(19,69.4) 



from which the allowed solutions are 

v = =>• n = 0, n= 2' n = 1 

z/ = ^ -, for all n except n = 2 (31) 

Our choices of parameters in preceding discussions are based on this argument and also 
constraints on non-minimal coupling from reference [19]. Especially, we have taken into 
account the following considerations: the solutions with v — are not interesting at all since 
they provide static cosmologies with non-evolving scale factor on the brane. In ordinary 
4D gravity with the action of the form F(R) = IqR! 1 and without any matter, a constant 
scale factor with a singular equation of state is achieved for n = 1. Similarly, in the normal 
DGP model, that is, when n = 1 and with no ordinary matter present, we cannot define 
the equation of state because the scale factor on the brane is constant and equation of state 
parameter will be singular [10,20]. In our F(R,<f>) setup with F(R,<j>) = \{l - £(j) 2 )l R n , 
as we have discussed previously, to have accelerated expansion we need to rstrict n to the 
interval 1.37 < n < 2. In this situation, there is no term with the first power of the Ricci 
scalar to obtain ordinary Einstein-like interactions between two astrophysical objects, such 
as the Sun and Mercury. To overcome this problem in our setup, we have treated the case 
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with n = 1 separately after equation (25) above. On the other hand, if we choose theories 



such as F(R,<j>) = \{1 - £<f) 2 ) 



where ( is a suitably chosen parameter, 



our model will contain the case with n = 1 and therefore general relativity. We will come 
back to this issue later with more details. 




Figure 1: Choosing rj = 2, £ = —0.166 and n = 1.4, we find z/ = 1.2 from equation (27). 
There are two intersection points accounting for positively accelerated expansion within F(R, 4>) = 
i(l - i<t) 2 )^R n scenario. 



4 Dynamics of Equation of state Parameter 



Because of explicit interaction between scalar field and curvature, the equation of state 
parameter for non-minimally coupled scalar field in this F(R, <fi) DGP-inspired setup is as 
follows 

p(tot) p(curv) _|_ p 



Jtot) 



p (curv) + ■ 



(32) 



which can be written as 



<j> 2 - V(j>) + 2 



2RH + R F"(R, 4>) + R 2 F"'(R, 0) - ±F(R, 0) + ±RF'(R, <j>) 



14,2 + vw _ 6 dF Hil+ _ 7 _a_ ( 1 F(Ri 0) _ i RF , (R 



dR 1 J 



(33) 
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Figure 2: Choosing rj = 4, £ = —0.166 and n = 1.4, we find z/ = 1.2 from equation (27). There is 
one intersection point accounting for equality of two sides of equation (27) and positively accelerated 
expansion with this choice of parameters. 




-20 J 



Figure 3: Choosing rj = 2, £ = +0.166 ( a positive non-minimal coupling) and n = 1.4, we find 
v = 1.2 from equation (28). There is one intersection point accounting for equality of two sides of 
equation (28) and positively accelerated expansion with this choice of parameters. 
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Figure 4: Choosing r/ = 4, £ = +0.166 and n = 1.4. In this case there is no intersection point 
between two sides of equation (27) corresponding to v = 1.2. It is interesting to note that even with 
other values of v > 1, there is no intersection point. Therefore, it is impossible to have accelerated 
expansion in this case. 

This is a general statement for equation of state in F(R, 0) setup. To proceed further, 
we have to specify the functional form of F(R,<p). Firstly we assume a conformal coupling 
between scalar field and induced i? n -gravity as F(R,(f>) = |(1 — £,(fi 2 )£oR n . The scalar 
field potential is chosen to be V(<f) = and then using natural ansatz a(t) w At" and 
4>(t) 7a Bt'L 1 , the following time-dependent equation of state parameter will be derived 

u(t) = | M S 2 t 2 " + '" 1 ^t 2 " - 2H (\2 V 2 - Gv) n {v.+ v)\ - X B"t 2 ' l +' ln + 2 + 2ie a (12v 2 - 6i/)"t 2 ( £B 
(12 v 2 - 6u) n t>"> + A^S"i 2 ^+ 2 " ) + m 2 3 t 2 " +4 - n+l "' ( 2(n - 1)(2„ + 2n - 1) - 3^(2k - 1)(- - 1) 



(34) 



-B 2 M 2 i^" , + 4 " + AB"t 2 *'+ 4 "+ 2 - 6£B 2 iiv(12v 2 - 6u)"t^+ 2n + ?, v (2v - - 1) + 6^n(n - i)i 2 ^+*"l+4r. 



The equation of state parameter in this scenario has a very complicated form. To find 
further intuition, this equation is plotted for several interesting cases in which follows. As 
these figures show, essentially crossing of phantom divide line, u — — 1, is supported in this 
F(R, 0) DGP-inspired model. However there are some points that should be stressed here. 
As figure 5 shows, the minimal case for a single scalar field has no phantom divide line 
crossing for some region of parameter space. On the other hand, this non-minimal F(R, 0) 
setup has no crossing of phantom divide line with negative n values. However, as table 2 
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and related arguments in preceding section show, accelerated expansion is possible even with 
negative n values. This is supported by other studies [21]. For n = — 1 we have a power-low 
acceleration on the brane without having to introduce dark energy. This result is consistent 
with the observational results similar to dark energy with the equation of state parameter 
— 1 < uj < — | which has no phantom divide line crossing. Figures 6 and 7 show the time 
evolution of equation of state parameter with some suitable values of model parameters. 

As we have stressed previously, the case with n — 1 needs further discussion. We have 
adopted the ansatz F(R,(f>) = |(1 — £0 2 )£o-R™ and our analysis has shown that n should 
be restricted to the interval 1.37 < n < 2 to have reliable cosmology in this setup. In 
this situation there is no term with the first power of the Ricci scalar to obtain ordinary 
Einstein-like interactions between two astrophysical objects, such as the Sun and Mercury. 
To overcome this problem in our DGP-inspired F(R, 0) gravity we have treated the case with 
n = 1 separately in subsection 3.1 and after equation (25). Form more general viewpoint, 
we can unify these arguments by choosing the following form of F(R, 0) 



F(R, 0) = hi - £0 2 ) \R - (1 - n)C 2 (5)1 , (35) 



2 



\ ( 2J 



where ( is a suitably chosen parameter [3,21]. The accelerated expansion and crossing of 
phantom divide line in this setup can be studied in the line of previous arguments in this 
paper. Fortunately this type of theories contain n = 1 and general relativity as a subset. 
For instance, the equation of state parameter in this setup takes the following form 

u,(t) = ji M 2 B 2 t~ 2M ~ 2 - AST™ +2iBADt-"- 2 ^Bi-"- 2 + S^St""" 2 + S"" 1 Xr,t~ + ^ + 4«S 2 AD^f- 2 "- 4 - HI_£^ M 2 B 2 



-2/j-2ti+1 m 2 A 



1 + n(D - 



2nd t ""(n - 1) (-2v + 3A* + In - 4) - -Dt 



J Vb 2 ^- 2 + AS"*-- - 6M.es 2 ADt- 2 ^ - ^I.^-Mn- l) 2 *- 2 " 1 
| 2 1 + n(D - l)t3»-l 

where A = Ylv 2 - 6u, C = C 2(1_n) and D — 1 — (1 — n)CA n - l t n . Figure 8 and 9 show 
the result of numerical analysis of this model. We see that these types of theories essentially 
account for crossing of phantom divide line. In this case as figure 9 shows, we have crossing 
even with negative values of n. Note that by choosing negative values of n, we can treat 
theories with F(R, 0) = 1(1 - £0 2 ) (r - ^g^) in this framework. The model proposed by 
Carroll et al [21], lies in this framework. For treating the issue of late time acceleration in 
this case, Friedmann and acceleration equation with this choice of F(R, 0) take the following 
forms 

3u 2 mlD(l-^B 2 r 2 ^) = ^Bh-^+2\B^r m+2 -l2^ABh- 2 ^- 2 +l2^CA n B 2 t- 2 ^ 2n (l-n) + 
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(mi - l)A + — CA n -\n - l) 2 m^r 2n+2 , (37) 

and 

-3i/(f-l)m3(l-5Dfl 2 r 2f ') 
6£B 2 AIJi~ 2 ' i ~ 2 (-£DA +3/jy 

where A, C and D are defined previously. Numerical analysis of these equations shows the 
possibility of intersection points and therefore accelerated expansion. For instance, figure 10 
shows the situation for a specific choice of the parameters. 

As another important point, the issue of stability of the self- accelerated solutions should 
be stressed here. The self-accelerating branch of the DGP model contains a ghost at the 
linearized level [22] . Since the ghost carries negative energy density, it leads to the instability 
of the spacetime. The presence of the ghost can be attributed to the infinite volume of 
the extra-dimension in DGP setup. When there are ghosts instabilities in self-accelerating 
branch, it is natural to ask what are the results of solutions decay. As a possible answer 
we can state that since the normal branch solutions are ghost-free, one can think that the 
self-accelerating solutions may decay into the normal branch solutions. In fact for a given 
brane tension, the Hubble parameter in the self-accelerating universe is larger than that 
of the normal branch solutions. Then it is possible to have nucleation of bubbles of the 
normal branch in the environment of the self-accelerating branch solution. This is similar 
to the false vacuum decay in de Sitter space. However, there are arguments against this 
kind of reasoning which suggest that the self-accelerating branch does not decay into the 
normal branch by forming normal branch bubbles ( see [22] for more details). It was also 
shown that the introduction of Gauss-Bonnet term in the bulk does not help to overcome 
this problem [23]. In fact, it is still unclear what is the end state of the ghost instability 
in self-accelerated branch of DGP inspired setups (for more details see [22]). On the other 
hand, non- minimal coupling of scalar field and induced gravity in our F(R, 0) setup provides 
a new degree of freedom which requires special fine tuning and this my provide a suitable 
basis to treat ghost instability. It seems that in our model this additional degree of freedom 
has the capability to provide the background for a more reliable solution to ghost instability 
due to wider parameter space. 

Finally, the phantom divide line crossing in conventional scalar field theory violates posi- 
tive energy theorems. In our setup, non-minimal coupling of scalar field and induced gravity 
has the capability to evade this problem. In fact, following Lue and Starkman [24], if we 
consider the FRW phase instead of the self-accelerating phase, and relax the presumption 
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- M 2 S 2 *- 2 ^ + AS^-^^+ 2 -6M^S 2 ADt- 2 ^- 2 -6-^ 1 ,n(n-l) 2 CA"- 1 t- 2 "- 2 + -M 2 S 2 ^ 2A '-3Ai3''t-^''+ 2 - 
2 D 2 

-A^S"- 2 i-'"'+ 2 ' i + 2 )-6M 2 «BAt- 2 ^- 2 + 6CM 2 CA"S 2 i- 2 ' i - 2 "(l-n) + -|( 2,,'A" '„,,. : " 

9A"~ 1 Cn(n - l) 2 i~ 2 ™ + 2 - -D, (38) 
4 



that the cosmological constant be zero (i.e., abandon the notion of completely replacing dark 
energy), then we can achieve uj < — 1 without violating the null-energy condition, without 
ghost instabilities and without a big rip ( see [24-27] for more details). 




n=1.37, nu=1.02, xi= -0.83 , eta =2 

n=1.39, nu- 1.13, xi=-0.83 ,eta -2 

n-1.37, nu- 1.02 , xi= -0.83, eta =4 

n-1.39, nu-1.13, xi- , eta =4 



Figure 5: Crossing of u = — 1 line in F(R,<fi) = ^(1 - £,<p 2 )£ R n DGP-inspired model. As the 
figure shows, with a single minimally coupled scalar field there are regions of parameter space with 
no crossing of phantom divide line. 

5 Summary and Conclusions 

Current positively accelerated expansion of the universe could be the result of a modification 
to the Einstein-Hilbert action. Also DGP braneworld scenario has the capability to interpret 
this late-time acceleration via leakage of gravity to extra dimension in its self-accelerating 
branch. On the other hand, there are several compelling reasons for inclusion of an explicit 
non-minimal coupling of scalar field and induced gravity on the brane sector of DGP action 
[14,16]. These arguments have led us to study cosmological implications of a DGP-inspired 
F(R, 0) gravity. Although this issue has been studied extensively for the minimal case, our 
setup provides a new approach to treat accelerated expansion within a general framework 
of modified scalar-tensor theories taking into account the role played by curvature correction 
and non-minimal coupling of scalar field and modified curvature simultaneously. We have 
studied late-time acceleration and possible crossing of phantom divide line in this setup. The 
condition for accelerated expansion and also equation governing on the dynamics of equation 
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xi= -0.123 , eta =2 




1 2 3 4 



t 



11=1.4, nu=1.2 — 


— n=1.5 , nu= 2 - - 


- - ii=1.6, nu=3.3 


n=-0.37, nu=1.005 







Figure 6: Crossing of uj = — 1 line. This figure shows that negative values of n cannot account for 
phantom divide line crossing in this non-minimal setup. 



xi- -0.123 




-4-| , , , , , , . , , 

1 2 3 4 5 



t 



ii=1.37, nu=1.02 , eta =2 


11=1.39, nu=1.13 ,eta=2 


n=1.37,eta=4 


11=1.39 , eta =4 



Figure 7: Crossing of u = — 1 line in F(R,<p) = ^(1 — £(j) 2 )l R n DGP-inspired model for different 
scalar field potentials and negative non-minimal coupling. 
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mu=0.9, xi= -0.166 




- n=-l , nu=2 , eta =2 

- n=1.4 , nu= 1.2 , eta =2 



■ n=- 1 , nu= 2 , eta =4 
11=1.4 , nu=2 ,eta =4 



Figure 8: Crossing of u = -1 line in F{R,<j>) = \{l - icp 2 )[R - (1 - n)C 2 (^) H 
model for different scalar field potentials and negative non-minimal coupling. 



DGP-inspired 



mu=0.9, xi= -0.166, eta -2 




- H--2 , nu=15/4 
n=1.5 , nu=2 



n=l .4 , nu= 1 .2 



- n=-l , nu= 2 



Figure 9: Crossing of u = -1 line in F(R,<f>) = \{l - £0 2 )[r - (1 - n)( 2 ^Y 
model for different n and negative non-minimal coupling. 



DGP-inspired 
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Figure 10: With rj = 2, £ = 0.166 and n = 1, we find i/ = 1.2 from equation (38). There is 
an intersection point between curves corresponding to two sides of this equation with v = 1.2 
accounting for positively accelerated expansion within F(R, </>) = 1(1 — £0 2 ) i? — (1 — n)£ 2 ( 
scenario. 



of state parameter are very complicated, so that we were forced to try a reliable and natural 
ansatz to find some intuition. We have shown that this setup accounts for accelerated ex- 
pansion with a suitable choice of parameters space or fine-tuning. Also, this setup accounts 
for crossing of phantom divide line in some ranges of parameters appeared in the model. In 
the minimal case, our model coincides with existing models of f(R) gravity. Especially the 
case with n = — 1 in the absence of non- minimal coupling exactly coincides with well-known 
result [21]. On the other hand, for a single scalar field with minimal coupling to gravity 
there are region of parameter space with no phantom divide line crossing. We have seen 
in the numerical calculations that with V((f>) = for both 7] = 2 and rj = 4 and with 
arbitrary £, there is no crossing of phantom divide line for negative values of n ( Figure 6 ). 
However positively accelerated expansion of the universe can be explained in this situation 
with negative n. On the other hand, if 1.37 < n < 2, for negative £ we have crossing of 
the phantom divide line ( see figure 7 ). As an important result in the analysis of late-time 
behavior, from equation (27) we find that only for positive values of n with 1.37 < n < 2 and 
for positive or negative values of £, there are intersection points and therefore accelerated ex- 
pansion for corresponding parameters choice. The case with n — 1 which describes ordinary 
general relativity has been treated separately. By adopting more general theories such as 
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F(R, 0) = i4> 2 ) R— (1— n)( 2 ^^j , the previous restriction on n is relaxed so that n — 1 
is a natural subset of the model. In summary, accelerated expansion can be explained by both 
positive and negative values of n with suitable fine-tuning of other parameters. On the other 
hand crossing of phantom divide line in F(R, 0) = |(1 — £<J) 2 )£ R n model can be explained 
only with positive values of n and fine-tuning of other parameters. The range of variation of 
n to have crossing of phantom divide line in this setup is 1.37 < n < 2. However in theories 



with f(r,<i>) = 



this crossing occurs for both negative and 
R-(l-n)C 2 (§) 



n 



R-(l-n)e(^) 
positive values of n. For theories of the type F(R, 0) = |(1 — £0 2 ) 
is not restricted to this interval and these theories account for phantom divide line crossing 
even for negative values of n as figure 9 shows. The issues of ghost instabilities and violation 
of positive energy theorems have been discussed also. Due to wider parameter space as a 
result of non-minimal coupling of scalar field and modified gravity, we hope this model can 
evade these problems. 
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